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The evolution of primordial black hole masses in the radiation-dominated era 



Paulo Sergio Custodio and J.E. Horvath 



We discuss in this work the behaviour of primordial black holes (PBHs) in the radiative era. 
Taking into account the Hawking evaporation and the absorption of energy we revisit the complete 
differential equation for the evolution of the mass of a PBH. We show that the mass can grow in this 
cosmological phase in a very slow fashion (even when considering the very high temperature of the 
radiation) if at all, and give a strong upper limit to the maximum accretion of mass. We evaluate 
relativistic effects due to the peculiar motion relative to the CMBR and show that the existence of 
relativistic black holes with very high mass absorption is highly unlikely. Finally we demonstrate 
that thermodynamical equilibrium between black holes and the cosmic radiation can not exist for 
finite times, and therefore initially non-evaporating PBHs must jump to the evaporating regime. 
This analysis supports the several efforts performed to look for signatures of evaporating holes. 



I. INTRODUCTION 



The search for PBHs is certainly one of the most interesting areas of high-energy astrophysics. While quite strong 
limits have been published [?] on the fractional mass in PBHs, it is possible that a fraction of short GRBs may 
correspond to evaporating holes having masses at formation of the order the Hawking mass M <~ 10 g that would 
be evaporating today [?]. The celebrated Hawking radiation would be driving black holes to their end, although 
several important details regarding this final phase are still being analyzed. However, Sivaram [?] has argued that 
these objects would not be observed at all, since that these PBHs would have absorbed a substantial amount of mass 
when the Universe was very young and hot, and therefore would not be evaporating today. This cautionary remark is 
important, since then short GRBs must be produced by other sources. Here we reanalyze the problem, and show that 
substantial mass gain did not occur at all (Section 2). The kinematic analysis of Section 3 indicates that a possible 
boost to enhance the mass absorption term is severely suppressed. The last Section is devoted to a demonstration 
of the impossibility of keeping thermodynamic equilibrium between the black holes and the radiation, thus justifying 
the previous discussion and the expected evolution of the mass. A slightly more general argument along these lines is 
given in the Appendix. 



II. EVAPORATION AND ABSORPTION OF MASS 



Hawking [?] , found in a seminal paper that black holes would loss their mass through a slow but irreversible quantum 
process, where virtual particles from the vacuum are converted in real ones created by the large gravitational energy. 
This quantum process gives rise to mass loss where the black hole had its temperature increased continuously. The 
mathematical analysis of this problem, showed that this emitted spectrum is thermal, with an associated temperature, 

he 3 10- 7 K 

bh 8irGk B M ~ (M/Mq) ' 1 ' 

Using the Stefan-Boltzmann law and the formula above, the evolution of the mass of the PBH is given by 

dM __A{M) 

dt ~ M* ' ( ' 

where A(M) counts the degrees of freedom of the emitted particles. Numerically, A(M) is given by 

A(M) ~ 10 24 g 3 s- 1 r ] {M) ; (3) 

with r](M) = 7.8 x D\{M) + 3 x D 1 / 2 (M), D 1 / 2 and D\ are the multiplicity of quantum numbers (spin, charge 
and color) of the emitted particles. For the standard model with three fermionic generations, we have D\ — 27 and 
D 1/2 = 90 at T ~ lOOGeU and A(M) < 7.8 x lO 26 ^" 1 (ignoring contributions from Higgs scalars and gravitons, 
sec [?] for details). 

We know that the cosmological environment was very hot and dense in the radiation era, therefore, we expect some 
classical absorption of the energy-matter from the surroundings of a given PBH. We adopt a Friedmann-Robertson- 
Walker model, with R(t) the scale factor describing the expansion (or contraction) of the universe and K = — 1, 0, +1 
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are the possible values of the curvature parameter, describing, respectively open models (K = —1 or K = 0) and 
closed models (K = 1). The dynamics of the FRW model is governed by the well known equations 



R\ 8ttG K 

r) =— q -r^ (4) 

R 4ttG, 

R = ~~ {g + 3P) 5 (5) 

without cosmological constant (which is not relevant for our discussion) . In the radiation-dominated era the universe 
was close to K = and the scale factor and the radiation temperature evolve as R(t) = RiX t 1 / 2 and T ra d(t) oc R(t) 
respectively. Numerically, the temperature is given by T ra d <~ lMeV(t/ls) and then the radiation density is given 
by 

Qrad{t) ~ 8.4 x 10 4 (t/lsr 2 3 cm- 3 . (6) 

We define as usual the cosmological horizon by the integral Rhor(t) = R(t) lim^o J t * T^py- The cosmological 
horizon exists when this limit converges for finite t. Taking the limit tj — > we have for the solution R(t) oc t 1 / 2 

R hor {t) = 2ct . (7) 

The mass contained in the horizon is given simply by the product Mh or {i) ~ ^-Rhor(t) 3 g(t). Using eqs.(6) and (7) 
we obtain 

M hor (t)~ 7.6 x I0 37 (t/ls)g. (8) 

This relation is very important to the problem, since we expect that PBHs would be formed with masses comparable 
to the horizon mass [?]. There is also a dependence with the adiabatic index of the background matter T = ^, an 
important feature for the collapse of perturbations associated to PBHs and generally overlooked. Since PBHs with 
initial masses larger than the horizon mass are excluded by causality, for a generic PBH with initial mass M was 
formed in the early universe its formation time must satisfy 

tf(M) > 3 x I(T 38 (M/g)s . (9) 

Therefore, from tf(M), we expect that the PBH mass will be below the horizon mass at subsequent times t > tf. 
In order to address the complete problem of the evolution of PBHs masses, we shall include a classical absorption 
term. This term must be included to the r.h.s. of cq.(2) and is of form 

(^) b =°9W)F rad (T) ; (10) 

where a g (M) — ^-r g 2 is the gravitational cross section for the capture of relativistic particles [?]. Since F ra d(T) = 
cQ ra d{T) is the radiative flux as seen in the rest frame of the black hole, the complete differential equation for the 
mass is given by 

dM A{M) 27^ 2 . 

-dT--^^ + ^^^ ad (T)M , (11) 

where we have assumed that M > Mpi anc k in order to avoid complicated corrections due to quantum gravity effects 
and neglected the back-reaction (see [?] and references therein for a quantum mechanical evaluation). 

The inclusion of the second term is a necessary ingredient for the thermodynamical description, since PBHs formed 
obeying the causality constraints with masses > Muaw before t <~ Is will experience an inward heat flux from the 
hotter surroundings. Consequently, these cold PBHs will gain mass-energy at the expense of the radiation. When 
the physical temperature of these PBHs equals to the radiation temperature, these objects will be in instantaneous 
equilibrium and M = at this instant. We can derive the mass for equilibrium to happen at any epoch solving M = 

[?] 

10 26 o 

M c (t) ~ ^ ; (12) 
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where To is the present cosmic temperature. Any PBH with mass greater than this value at t must be accreting 
(M(M > M c ) > 0) and not evaporating, since the Hawking evaporation is negligible in this case. 
To be quantitative we may consider the ratio of the absorption to evaporation terms to obtain 

Mabs /M X 1 



« M) =iC=teJ ; (13> 

which shows that the Hawking term quickly becomes negligible for masses above M c . 

A first look at eq.(13) could convince us to accept the conclusions of Sivaram [?], namely that these supercritical 
(M > M c ) PBHs would be growing very fast in the early universe. To show that this would be misleading we shall 
solve the eq.(ll) for initially supercritical PBHs, (that is, ignoring the first term) in the thermal radiation bath 
evolving as stated in eq.(6). Then, the differential equation to solve is 

dM 27ttG 2 , . 

= —^Qrad(t)M 2 . (14) 

The integration of this equation is immediate and we obtain 

1 1 



— + 1.2xl0" 39 kT 1 ; (15) 

M(t) Mi \t Uj y v ' 

where Mi is the initial mass (at formation) and ti is the formation time (measured in seconds) for this object. This 
solution is valid if the initial mass satisfies M c (ti) < Mi < Mh or (U), as stated above. 
We may invert the analytical solution cq.(15) above, casting it in the form 

with e(M,,( i )= 1.2 xl(T ss g!jM 

Since e(Mj, ti) << 1 for PBHs satisfying M, << Mhor(U)', and even in the most favourable circumstances e(Mj, ti) ~ 
Mhor{ti), we can write eq.(16) expanding it in a Taylor series and find that the growth of these objects is M <~ 
Mi[l + e(Mj, U)]. Note that the mass gain is larger for larger masses, and also increases for early formation times for 
a given mass. In fractional terms 

SM M — Mi 

w = ^r=<M i ,t i ). (ir) 

Note that the fractional mass gain saturates with time for t >> ti, as expected from the time dependence of the 
radiation density Q ra d(t) oc t~ 2 . Even at t = ti the universe was not dense enough to drive a very strong growing. We 
can show it in the following way. Suppose that the hot environment can indeed support a relativistic PBH mass gain 
of the gravitational radius r' g ~ c. Let us show that the actual mass accretion was many orders of magnitude below 

this value, i.e. M « ^~ 2 x 10 38 gs _1 . For this purpose, assume that initially, we had M — 27 *p 2 QradM 2 ~ 
Solving it, we find the minimal density of the radiation capable to drive the relativistic growing for that mass M 

Qmin - 5.7 x 10 51 gem- 3 (M/M Haw )~ 2 . (18) 



If we compare Qmin with eq.(6) we obtain that these PBH would have to be formed before t re i(M) 



x 



3.8 x 10 {M/MHaw)s. However, according to eq.(9) causal PBH formation is allowed after tf > 3 

10 -23 [M/M Haw) s. Thus, most of the time the actual growing has to be much smaller than the maximal theoretical 

3 

value M m ax = §q , which needs very extreme conditions. 

A brief relativistic gain may nevertheless be allowed, since tf and t re i are very close. This means that a PBH with 
horizon mass, formed exactly at tf will gain some mass, but the maximal fractional gain is limited to be < 0.04. 
Therefore, if we consider the maximal gain, a PBH formed with Mi ~ 10 15 g would finish its evolution with final mass 
(at the end of the radiation era) greater than the Hawking mass, but with a small gain. However, that is enough to 
delay its evaporation, since the evaporation timescale is given by 

tevap(M) ~ t (M/M Haw ) 3 ■ (19) 

therefore the mass M ~ M 4 + 0.04Mi then we obtain t evap (M) ~ (1.04) 3 t - 
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The situation is quite different for PBHs with initial mass of order Mj smaller than Mn aw . In this case, even the 
maximal gain can not drive these PBHs above the Hawking mass at all. Therefore, these objects will show very small 
modifications of their behaviour. Note that the maximum gain is easily obtained for Mj <~ M^ or . Then substituting 
we have for that absolute maximal gain 



^-rnax 

{Mi -> M hor ) ~ 1.17 x 10" 39 M; 



1 



3 x l(T 38 M f 



< 0.04 . 



(20) 



In order to understand the full behaviour of the PBHs masses it is useful to consider the evolution of the critical 
mass. The latter may be expressed in the form 



M c (t) ~ 10M Haw (t/ls) 



1/2 



(21) 



If we compare the behaviour of a PBH with initial mass above the critical mass (but obviously below the horizon 
mass) and the growth of the critical mass, we observe that this object will cross the critical mass some time after its 
formation. While this PBH stays above the critical mass, it will not evaporate, but its final mass (that is, before the 
crossing time) will be given by eq.(19). The crossing time t cross is defined by 



M(t) = Mi[l + e(Mi, U)] ~ 10M Haw {t cross /ls) 



1/2 . 



with solution 



,(Mj) ~ 10~ 2 s(Mi/M Haw y[l + e(Mi,U)Y 



(22) 



(23) 



This explains why these objects did not gain much mass: the time interval in which these objects stay above the 



critical mass is very small if compared to cosmological times, or to the timescale for evaporation t e 



t M 3 . 



Any supercritical PBH will be driven by the second term in eq.(14) and any subcritical PBH will be driven by the 
first term, the Hawking evaporation. In the section 3 we will show that the objects crossing the critical mass curve 
can not stay there, and must drift to the evaporation region of Fig.l in which the evaporation quickly dominates. 



III. MASS GAIN INDUCED BY HIGH PROPER MOTIONS 



Since we expect that at least for large Lorentz factors the absorption and evaporation processes must be modified, 
we want to consider corrections due to the proper motion of the PBH relative to the background radiation. This effect 
is described by the following expression 

lrad{V,V) - —g— , (24) 

7(1 - n.P) 

where we will denote T ra d(9,v = 0) simply as T ra a, the temperature of the background radiation in the rest frame 
of the PBH when this PBH has v = 0. Then, T rac i(0,v) for arbitrary v, denotes the angular distribution of the 
background radiation seen from the rest frame of this PBH when it moves with peculiar velocity given by v. The 
angle 6 is defined in the scalar product n.(3 — v[3cos(9). n is the vector defining the sight of view for an arbitrary 
observer located at the Schwarzschild radius, 7 = ^ and (3 — (v/c). 

The radiation density has an angular distribution as seen from a stationary observer located at the black hole 
horizon. The sign in the denominator of eq.(24) is such that we have a blucshift in the direction of the motion and a 
rcdshift in the opposite. The absorbed radiation is an angular average and it is given by [?] 

< Qrad(T r ad) > = J dQ.Q rad {T rad {6 , </>)) . (25) 

This integral is performed to give 

4 7 2 - Qrad(Trad(t)) ] (26) 
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where the time dependence in T ra( i(t) arises from the cosmological expansion, i.e. T ra d oc R . 

Inserting this result in the eq.(ll), we arrive at the differential equation for the mass of a black hole moving with 
velocity v, as seen by an asymptotic observer. Evidently g ra d in eq.(26) must be the same quantity given by eq.(6). 
Now, we seek a 7 capable to drive the rate of mass accretion to the maximal rate M 1 ~ Scaling the black hole 
mass to Hawking mass units 



1 



7 x 10 5 



4 7 2 - 

1) {M /Ad Haw)' 

This is a second degree equation for 7 with solution 



■{Qrad{t)/gcm~ 



7(t) 



8/I 2 Qradit) 



i + v /i + i|M 4 g rad (t) 2 



(27) 



(28) 



where a = 7 x 10 53 , [i = M/Muaw and g ra d(t) is expressed in gem 3 . Substituting g ra d(t) = g(tf)(t/tf) 2 , we can 
rewrite eq.(28) as 



7(*) = 



<t/t f y 



8/j, 2 g r ad(tf) 



16 

1 + \l 1 + T2^Qrad(tff(t/t f y 4 



(29) 



since g ra d{t) ~ 8.4 x 10 4 (i/ls) , and our solution holds only for t > t/(M) we obtain the following inequality 

10 51 



Qrad{t > tf) < 4.4 X ^ 2 



(30) 



Therefore, the prefactor in eq.(28) has to be greater than <~ 18, for any t> tf. We may rewrite the required Lorentz 
factor as 



7(t) > 18 



16 2 

1 + -^H 4 grad(t) 



Applying the same arguments to second term, we can finally write 



7 (t) = m/tfY 



l + ^/l + 6x 10- 4 (^/)" 



(31) 



(32) 



For asymptotic times t > tf, we have j(t) — > 36(t/tf) . 
Therefore, we conclude that even at the beginning of its life (t 



t/) and when the PBH was formed with M ~ M; 



/tor ) 

a high value for the proper velocity v would be needed to drive the PBH into the fast growing mode (differing from 
c in one part per thousand) . The additional terms in the square brackets are positive and enhance the value of the 
minimum 7 in eq.(31). The physical reason is very simple: from the initial time formation, the radiation density 
decreases with time, and consequently we will need to impart very high values for the proper velocity to drive an 
object above the critical mass. This particular situation was certainly unlikely in the early universe, and we do not 
expect to have PBHs with substantial accretion of mass. It is worth to remark that the ingredient that allowed 
us to perform these numerical evaluations and conclusions is the concept of critical mass (see also [?] for a related 
discussion). 



IV. IMPOSSIBILITY OF THERMODYNAMIC AL EQUILIBRIUM BETWEEN PBHS AND THE 

BACKGROUND RADIATION 

In this section we address the problem whether black holes and the surrounding radiation can be in thermodynamical 
equilibrium. This would imply that the PBH mass to be constant. Page [?] analyzed this problem in a gedanken 
experiment with one black hole contained in a finite box plus thermal radiation at a constant volume. It is easy to 
understand that this black hole within this box will settle down in stationary equilibrium, with its mass remaining 
constant in time. In this case, the total energy of the radiation + black hole system is 

E = AT 4 V box + M bh . (33) 



6 



Since that this system is closed, we have for variations in mass and temperature 



SM = -4AV box T 3 ST . (34) 

If the black hole loses mass due to Hawking effect SM < 0, then the radiation temperature must raise ST ra d > and 
drive the black hole into equilibrium. If the black hole absorbs too much mass the radiation temperature diminishes 
substantially and additional accretion will not be allowed. We see immediatly that the situation is more complicated 
if the box has variable volume. It is immediate to conclude that microscopic variations in the critical mass must obey 

<5M c cx (T rad y 5 ^- . (35) 

Vbox 

For very high temperatures, the variations in the critical mass are very small, as expected. 

We shall now analyze the more complex and realistic situation, consisting in black holes and radiation within a 
cosmological model. For this, we emphasize that the critical mass is a cosmological quantity, defined by the radiation 
temperature. Numerically, we have 

Therefore, since the temperature decreases with cosmological time, the critical mass grows monotonically. However, 
we note that the critical mass allows an instantaneous equilibrium, since that M(M — M c ) = 0. As the universe 
expands a black hole initially in equilibrium with the background radiation must deviate from it because of the 
different temporal behaviours. 

To prove this statement, let us evaluate the variation of the area of the black hole. We know that R g = ^r^, then 
the horizon area is given by A = 4irR g 2 , and therefore we have 

3?7T 

dA = °—G 2 MdM . (37) 

c 4 

The sign of the mass variation is related to the energy that crossed the event horizon as dE = c 2 dM. Then an 
energy flux through the event horizon induces a variation in the area given by 

^ = C * (W) 
dA 32nG 2 M ' k ' 

Then, the time derivative of the energy flux per unit area will be given by 

/ d 2 E \ c 6 dM 



\dtdAJ 



32ttG 2 M 2 dt 



(39) 



Substituting eq.(ll) into eq.(39), we obtain 

r- 6 r a i 

(40) 



( d 2 E\ c 6 
{ dtdA ) ~ 32ttG 2 



A 27nG 2 

Jfi ^3 Qrad{t) 



We will denote eq.(40) as <r(M, t) for simplicity. We see that the energy is zero when M = M c , then a(M = M c , t) = 
for any t. However, by evaluating the time derivative at M = M c and we will prove that the first derivative at the 
critical mass is not zero, indicating an inevitable deviation from the initial equilibrium. 

Note that u(M, t) is positive-defined for hot black holes (subcritical black holes) and negative in the opposite case. 
Now, the time derivative of a(M,t) at the critical mass is given by 

27 

&(M = M C7 t)^~—c 3 g rad (t)>0; (41) 

and it is positive-defined since that the radiation density decreases with time. If we integrate eq.(40) with the 
initial condition M = M c (ti) until t > U and given that the radiation diluted during this interval due to expansion, 
we obtain the result 



ra(M c ,U) c 6 rt 

JU) da = -^G 2 J ti dt 



AA dM 27nG 2 . . 



(42) 
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and the right side can be rearranged in the form — Jm 7$ ~ fi c3 fo™^) ^Qrad- Then we arrive at 

Ac 6 



a(M,t) = a(M = M C ,U) 



327rG 2 



1 



32 
1 



' Brad(ti) 

27 3 A 

32 r 



(43) 



The first term a(M — M c ,ti) — by definition, and the last term is negative as we already know. If we consider 
that M — M c + SM, the term in brackets can be evaluated to give 



1 1 


1 


M A M c 4 _ 


" M c 4 



-1 



1 



1 + F{M C ,SM C ) 



(44) 



where F(M C , SM C ) = {^-f + 6(j^f + 4(^-) 3 + 4(^£). Therefore, the term in brackets is always negative, for 
arbitrary variations in SM, after the PBH had crossed the critical mass at t%. Note that F is negative because the 
terms containing higher powers of ^ are positive and they dominate the other terms (which are negative, since SM 
is arbitrary). For large variations of SM of any sign, we have F < 0, then we deduce that a(M,t) is positive, and 
thus conclude that this PBH migrates to the evaporating regime. 

Now, even if the variation SM is assumed to be small and negative, we could imagine that F(M C , SM) may change 
its signal and become positive. We shall prove that this is never the case. If the net variation in M is small, we can 
ignore the first term in the r.h.s. of eq.(40) and integrate to give 



a(M, U + At) 



0rad(ti)St ; 



(45) 



where St is a small time interval satisfying St <C 2771 "G 2 c 3 (M)Ag rad with g(M) 



AA 2 /M 7 - ABg rad 



In order to show that a PBH in an initial equilibrium M = M c deviates from it is enough to show that the first 
derivative & is non-zero at M = M c and is always positive for arbitrary variations SM. We proceed by steps. 

Since that we are interested in the algebraic sign of a(M, t) only, we change variables to absorb the coefficient 3 ^ G 
in eq.(40) and rewrite a as 



a(M,t) 



A 

M~* 



F>Qrad(t) 



(46) 



To prove that &(M = M c , t) > it is enough to see that its derivative is just &(M, t) — — M — B Q ra d(t) ■ However, 
the first term is null at M = M c and the second term is positive because we assumed that the universe expands, 
then &(M,t) > beyond the critical mass. Suppose that the variation in M is such that we write M — M c + SM 
and the "final" mass is below the critical mass by a significant amount, say M . Therefore, the energy flux conserves 
its original sign from the initial equilibrium. Now, let us prove that the sign of a(M,t) does not change if the small 

Bg ra dM 2 I — Bg rad {t). If we substitute 



A 

M 5 



A 
17P 



variation SM is negative. For this purpose, we write &(M, t) 
M = M c + SM taking into account that this change is very small, we can expand this expression and arrive at 

A 



&(M, t) 



A 

"M5 



M c 



+ 2M c B Qrad (t) 



SM > . 



(47) 



Therefore, the rate of change of a(M,t) drives the hole to the evaporation. To see that this condition of eq.(47) 
is enough we may calculate §f (M c + SM) > 0, thus implying that da(M c + SM) > or a(M c + SM) - a(M c ) > 0. 
This in turn means that a(M c + SM) > 0, since a(M c ) = and a(M) > happen in the evaporation region, or when 
the PBH losses much more energy than it is able to absorb from the environment. In the opposite case, if the mass 
variation is large but negative, we have a(M = M c + SM) = — Bg rad (t) <~ > 0, and the first term dominates 
because £ << 1 in this case. This completes the analysis for all negative variations SM around the critical mass. 

Now, we have in mind the general behaviour of the PBH mass evolution as the universe expands. From the eq.(46) 
we deduce that 



4A ■ 

a(M,t) = -—M - Bg rad (t) . 



(48) 



For subcritical black holes, M < and we know that g ra d(t) < due to cosmological expansion. Then, for 
subcritical PBHs is immediate that a > 0. 
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In the opposite case of supercritical black holes, the first term is —j^M < 0. Then, the algebraic signal of &(M, t) 
is determined by the difference —Bg ra d{t) — j^-M. However, we found in section 1 that the relative mass gain is 
strictly limited to a few percent. This means that the second term is negligible and the algebraic signal of &(M,t) 
remains positive. 

The algebraic sign of a(M,t) is always positive in "natural" circumstances. Some paths where &(M,t) < are not 
a priori excluded, but these paths correspond to very special conditions where g ra d > or the negative variation in 
Qrad{t) is very small. For realistic models of the radiative era, we know that g ra d(t) oc t~ 2 , then, this variation is not 
small compared to PBH mass gain. All those paths where &(M,t) < are entropically disfavoured. 

In the general case when the universe expands, the natural evolution for a PBH in the radiative era is to follow the 
type of (almost) horizontal path drawn in the Figure 1, before crossing the critical mass curve. There is a small mass 
gain along these paths, as we explained previously. The paths are thermodynamically favourable and terminate at 
the instantaneous equilibrium point M c (t — t c ) at the crossing time t c (Mi) (marked with a cross). Now, since that 
a(M,t) > at the critical mass, the deviation of equilibrium is unavoidable. The PBHs follow the type downward 
path afterwards, where a(M 7 t) > and its derivative are positive as already discussed. In other words, the Hawking 
evaporation will eventually dominate the classical absorption. 

In the general case we should analyze the full behaviour of 

27c 3 

- — Ag rad {t) ; (49) 



where M(t) is a solution of M = —jj?+Bg ra d(t)M 2 . If the initial conditions are given by Mi > M c (ti), g ra d(t) oc t~ 2 
(and it is implicit that g p bh(t) « g ra d(t) as before); and substituting the appropriate solution for M{t) into eq.(49) 
above, we deduce afte taking t » t c (Mi) that <j(M(t >> t c ),t) > 0. Therefore such a black hole crossed the critical 
mass curve and is now below the critical mass at t. The hole is then hotter than the environment and must evaporate. 

All these thermodynamical considerations show quite generally that an explosive growth of PBHs in the radiation 
era is forbidden. This conclusion in unavoidable when the problem is analyzed using carefully the critical mass. The 
validity of Sivaram's solution would also lead to the conclusion that an initially sub-horizon PBH would have to cross 
the horizon in a short time given by At ~ ^sZio^g''^ s ' This kind of behaviour was already ruled out by a set 
of calculations [?] using full General Relativity. Independently of the present results and former work, we are able 
to prove that this strong growth is impossible if Bekenstein's limit [?] holds (P.S. Custodio and J.E. Horvath, in 
preparation) . 

V. CONCLUSIONS 

We have analyzed in this work some features of the evolution of PBH mass in the radiation-dominated era. We 
found that within FRW models with null curvature there is no room for an explosive growth of their masses, which 
would need relativistic speeds of the gravitational radius. This result is in line with the previous analysis of Refs. [?] 
and [?]. Indeed the very modest value of the maximum fractional gain (0.04) argues against the possibility discussed 
in Ref. [?]. This means that essentially all the methods applied to limit their abundance make sense since evidence 
for these objects could lead to very valuable insights. We have also addressed in some detail two ways to preclude 
evaporation: a large proper motion altering the absorption-emission properties of the semiclassical PBHs and the 
possibility that they could be "dragged" and stay on the critical mass curve. In the first case we have given an 
ultrarrelativistic motion of the holes at their very birth (7 > 36) ; and demonstrated that the second situation is one 
of unstable equilibrium, and therefore that PBHs do cross the critical mass curve as discussed in [?] to evaporate 
(Fig.l). Other forms to avoid the onset of the evaporation might in principle be possible, but are not known to us. 
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VII. APPENDIX: A GENERAL ARGUMENT AGAINST A LINEAR GROWTH 



We shall now discuss further consequences of a linear growth of the mass as derived in Ref. [?] . If we assume that it 
is correct and substitute M(t) — Mi{t/U back into M(t) = B g ra d(t) M 2 (t) and impose that g ra d{t) = Q ra d(ti)(ti/t) 2 , 
we obtain for the r.h.s. 

—j£ = Bg rad (ti) = constant ; (50) 

for arbitrary times. Using again the same expression M(t) — Miitjti) and rearranging yields 

^ = BMf Qrad {U) ; (51) 

which shows that a linear solution would require a definite value of the radiation density at the formation time. 
Since initial black hole masses and the environment radiation density must be independent, we conclude that a linear 
growth is unphysical. 

We can generalize the above fine tuning constraint for arbitrary conditions provided reasonable physics still holds. 
Suppose that we do not know the behaviour of the radiation density at all. Then, it would be reasonable to guess for 
it some function of the form g ra d(t) — Qrad(U)F(t/ti). 

If we substitute this guess into eq.(17), we obtain 

/u 
du'F(u') ; (52) 

where u = (t/U). A power-law ansatz F(u) = u m for some m gives 

(1 - l/u) = ^-^[BM t g rad (U)}(u {m+1) - 1) . (53) 

The solution for m = — 2 displays the same unphysical fine tuning as before, since that it requires BMiQ ra d{ti) = 1. 

On general grounds there is no physical solution for m different from —2 which does not require some fine tuning 
between the initial mass Mj and the radiation density at this moment. Therefore, we conclude that the solution 
M(t) cat is ruled out. 
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VIII. FIGURE CAPTIONS 



Figure 1. Evolution of PBHs in the absorption and evaporation regimes 

The type of paths of a generic PBH are sketched in the PBH mass vs. cosmological time plane. The PBH is formed 
at the point marked with the star and undergoes a very small mass gain as long as its mass is larger than the critical 
mass (solid line), thus following an almost perfect horizontal path (long dashed) until the crossing time at the point 
marked with the cross. Since, as discussed in Section 3, the equilibrium is instantaneous, the PBH must drift to the 
evaporation region where the Hawking radiation will cause a decrease of the PBH mass and thus an evolution along 
the downwards bending path as indicated. Note that the mass as a function of time has a maximum on the critical 
mass curve; given that a has the inverse sign than M. The path claimed in Rcf. 3 is depicted with a dotted line and 
would need unlikely conditions as discussed in the text. 
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